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Abstract. Using the Fermi Golden Rule analysis developed in 
|CMj , we prove asymptotic stability of asymmetric nonlinear bound 
states bifurcating from linear bound states for a quintic nonlinear 
Schrodinger operator with symmetric potential. This goes in the 
direction of proving that the approximate periodic solutions of the 
NLS in [MWl do not persist for the fuU NLS. 



1. Introduction 
We consider the quintic nonlinear Schrodinger equation (NLS): 

(1.1) iut = -dlu + Vu- \u\%, m(0, x) = Uo{x), (t, x) G M x M. 

We assume the following hypotheses. 

(HI) The discrete spectrum of —d^ + V is ad{—dl + V) = {—Uq, —oji} 
such that Uq—Ui is as small as we want and is not a resonance, 
i.e. if M G L°°(]R) satisfies u" = Vu, we have u = 0. 

(H2) The potential V{x) is even: V{x) = V{—x) for all x and V{x) 
is real valued. 

(H3) The potential V is smooth and (x)"\/(^')(x) G L°°(M), for any n 

and any k, where (x) = \/T+~x^. 
(H4) Let ipj be real valued generators of ker(— (9^ + V + uj) with 

||'?/'j||i2 = 1. For {f,g) = f^f{x)g{x)dx we assume that for a 

fixed ao > 

(1.2) 5{r„^l)-{^ll)>ao>0, 

(1.3) 5{^'„^iy - {ro,l){^l^t) X^o. 

(H5) We assume that Fermi golden rule hypothesis (H5) stated in ^ 

Remark 1.1. The very strong regularity and decay hypotheses on the 
potential V{x) in (H3) are certainly unnecessary. See for example the 
dispersive estimates for for e^*^"^^"'"^-' with V G L-^'^(M) in [GScj, or 
the case with delta functions in |DMWj . Nonetheless, we do not try 
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to prove systematically the estimates stated later in ^ for these less 
regular potentials. 

We refer to the Appendix in |KKSW] about the existence of double 
well potentials satisfying (H1)-(H5), as well as to the brief computa- 
tional discussion in Appendix |B} Specifically, if one starts with an even 
potential Vo(x) such that —d"^ + Vq admits exactly one eigenvalue — fi, 
then setting 

V{x) = Vl{x) = Vo{x -L) + Vo{x + L) 

for L ^ 1 yields potentials with two eigenvalues, both very close to 
—Vt. Furthermore, for i~p{x) a normalized ground state for —d^ + Vq, 
then 

ijj{x) ^ 2-5 ((^(x - L) + {-y^{x + L)). 

Then, 

because they are all about 2T2||y9|||6. Originally, results of this nature 
appeared in the work of E. Harrell [H] . 



The equations (1.1) with V{x) = Vl{x) ioi L ^ 1 are the focus 



of recent research |KKSWl IMW] because of the rich patterns detected 



at small energies. The references |KKSW[ IMW] both treat (O) with 
replaced by a cubic nonlinearity (which for our purposes is a 
difficult problem due to the subcritical nature of the nonlinearity). The 
result in |KKSWj proves the existence of a family of nonlinear ground 
states e^^'^(j)^{x) of (1.1) that bifurcate out of the linear ground state 



e °ipo{x). For to > coo close to ojq the e (f)i^{x) are orbitally stable 



see §2.3 and even in x. At some critical u = u* the ground states 



bifurcate for u > u* in two families, one formed by even functions, 
which are unstable, and the other formed by non symmetric functions, 
which are stable. The arguments in |KKSWj are quite general, but 
here we double check that this behavior continues to hold also in the 



case of the quintic NLS (1.1) 



In |MWj , the existence of more complex long time patterns is ana- 
lyzed by studying the dynamics of a simplified finite dimensional sys- 
tem, obtained by selecting a finite number of variables of the NLS in 
an appropriate system of coordinates. Over long times it is shown 
to be a good approximation of the full NLS. In particular, this finite 
dimensional approximation of the NLS admits a larger class of time 
periodic solutions than just the standing waves. The question then 
becomes whether or not these new periodic solutions persist also for 
the full NLS equation. In |MWj it is conjectured they do not persist. 
In this paper we do not address the solutions considered in jMWj . but 
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nonetheless for an easier problem we provide the mechanism by which 
the full NLS disrupts periodic solutions of a simplified system similar 
to that in |MW] . In Appendix |A] however, we will present evidence 
that indeed similar dynamical solutions exist that would collapse via 
the asymptotic stability analysis presented here to a nonlinear bound 
state asymptotically. Such dynamics have been abstractly studied in 
jKKP] as well. 

We recall that |MWj simplifies the NLS by first choosing as system 
of coordinates the spectral decomposition of —d^ + V , and by then set- 
ting equal to the continuous components. Here, we consider instead 
a natural representation of the portion of if^(M) near the surface of 
asymmetric ground states. There are then natural finite dimensional 
approximations of the NLS admitting periodic solutions. They are as 
legitimate approximate solutions of the NLS as those in |MWj . although 
here we do not try to check as in |MWj if they are good approximate 
solutions. Our solutions are relatively easy because they live arbitrarily 
close to the surface of asymmetric ground states. 



When the full NLS (1.1) is turned on, these approximate periodic 
solutions do not persist because the ground states are asymptotically 
stable. Hence the periodic solutions of the simplified system, now split 
into a part converging in H^{M.) to the orbit of a ground state, and 



another part which scatters like free radiation (see Theorem 1.2). This 
part of the paper fits easily in the framework of the literature of asymp- 
totic stability of ground states initiated in jgWTl [5W2l [BPTl [BP2] . We 
recall that the most general results are in |Culj . which contains a quite 
general proof of the so called Fermi golden rule. In the present paper 
though, we treat a quite special situation, due to the hypothesis that 
ad{—d1 + V) consists of just two eigenvalues, and so it is enough to use 
the simpler framework of |CMt ICu3] (we recall that |Cu3j is a revision 
and a simplification of |Cu2] . which contains various mistakes). To ad- 
dress the solutions in |MWj one can probably proceed similarly, though 
the complexity of the dynamical systems studied and the cubic non- 
linearity makes the analysis rather challenging. The difficulty though, 
is that the solutions in |MWj . while of arbitrarily small energy, might 
nonetheless not be sufficiently close to ground states. The issues then 
seem in some sense more "global", and closer in spirit to the problems 
addressed in jTYtlSWH] . 

Following |KKSW] . we consider nonlinear ground states of the form 

e'"*(coV'o + CiVi + ^?(a;,i))- 

Applying the machinery of |KKSWj . we prove that there is an uj* > Uq, 
with u* — Uq ^ Uq — Ui, such that for u G {coq, u*] there is a uniquely 
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defined map a; — )■ 0^ in C^{{uo,u*], H'^'''^(M.)) for any {k,s), wliere 

(1.4) WuWh'^.s := \\{xy{l-dl)^^u\\L2. 

For u > u* there is a bifurcation, with a branch of even ground states, 
and a branch of asymmetric ground states. We focus on the latter ones. 
We then prove the following 

Theorem 1.2. There is a 6o > such that for any u G {u*,uj* + 6o) 
there exist an eo > and a C > such that if 

\\uo - e^'^^cullffi < e < eo, 
there exist uj± G {lo*,lo* + 6o), 6 e C^(]R;]R) and h± G with 

\\h±\\m + \<-^± — < Ce 

such that 

(1.5) lim \\u{t, ■) - e'^^'Uu± - e''^'h±\\Hi = 0. 

It is possible to write 

u{t, x) = e'^^^^(f)^(t) + A{t, x) + u{t, x) 

with \A{t,x)\ < CN(t){x)^^ for any N, with lim\t\^oo C n (t) = 0, with 
\imt^±ooUj(t) = u±, and such that the following Strichartz estimates 
are satisfied: 

(-'-•6) ll'"llLt°°{R,Hi(R))nLf(R,Wi'"\R))nLf(IR,Lg°(R)) — 

Once the necessary spectral hypotheses in \CM\ ICu3] are proved in 
Section |3| Theorem 1.2 is a direct consequence of |CMt ICu3] . Nonethe- 
less we give a sketch of the main steps in the proof. In particular we 
review in Section |4]the material on dispersion of linear operators needed 
in the proof. Here we recall that the absence of the endpoint Strichartz 
estimate on M requires some surrogates. The surrogates were found by 
Mizumachi [Mj- However it turns out that ^ can be substantially 
simplified, and that the smoothing estimates contained in [M], while 
interesting per se, are not necessary in the proof of the main result 
in [M]. In fact the classical smoothing estimates introduced by Kato 
in [K] are sufficient. This is discussed in |CTt ICu3] and is reviewed 
in Section |4j See also the recent results of |DMWj to allow singular 
potentials in our analysis with restrictions to /c < 1. 



instability of approximate periodic solutions for nls 5 

2. Ground states for u starting at uq 

2.1. Ground states for uj starting at loq. As in |KKSW] we consider 
ground states of tlie form 

with po and pi in M and for r]{ipo, ipi,uj) a real valued function belonging 
to iJ^'''(]R, M) for any {k,s) with {r],ipj) = for j = 0,1. We are 
looking for the simplest asymmetric ground states possible and not for 
all possible nonlinear ground states branching out of the linear ground 
state. 

We denote by Pc the projection to the continuous spectral component 
of —dl + V. Hence we look at the system 

- uopo + ujpo - {ipo, (poV'o + pli'i + vf) = , 
(2.1) - uipi + upi - (^1, (poV'o + Pii'i + vf) = , 

{-dl + V + u)r] = Pcipo^o + pi^i + V? ■ 

By an elementary application of the implicit function theorem one ob- 
tains the following 



Lemma 2.1. For po and pi sufficiently small, the third equation (2.1) 
admits a unique solution rj = ri{pQ,pi,u) which depends smoothly in 
(po,pi, w) with values in H''''^{R,M.) for any {k, s) and can be expressed 
as 

4 

(2.2) r]{po,pi,uj) = ^pI~^ p{rij{po,uj) + plE{po, pi,uj), 

j=0 

where rij{pQ,u){x) = {—iyrij{po,uj){—x) for all x. 



The proof of Lemma [2TT] is a standard application of the impicit func- 
tion theorem and a resolvent identity. Similar expansions are proven 
in Propositions 4.1 and 4.3 in |KKSWj . 

Lemma 2.2. There is a fixed number Eq > such that for po G [—£o, Eq] 
admits a unique function u = uj{pq), with uj{0) = coq, such that 

(Po,0,r7(po,0,a;(po)) 

is a solution of system (2.1 ) in C°°{{—eo, Eq), H^''^{R, M)) for any {k, s). 



Proof. For pi = one can see that the third term in the Ihs of the 



second equation in (2.1) is 0, because it is /jg'?/'i(po'?/'o + Po'7o)^'^^) which 



vanishes since the integrand is an odd function. So the second equation 
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in (2.1) is trivial. Substituting r7(po, 0, w) = plrio{po,u) in the first, and 
factoring out a common factor po, we get 

-uo + u- po{ipo, (ipo + pli'i + ptvoipo, (^)f) = 0. 

By the imphcit function theorem we can solve with respect to u getting 

(2.3) u = u{po) =Uo + pti^l 1) + 0{pI). 

□ 

Lemma 2.3. There is a number pi e (0,£:o); with p^ ^ {ujq — uji)^ such 
that at Pq = the function of Lemma 2.2\ satisfies also the equation 
(2.4) 



Pi 



lpi=0 



Proof. Equation (2.4) is, for rjj = rij{po,u{po)), equivalent to 

d 

-ui+ u{po) - ^|p,=o (^1' (/'o^o + Pi^i + vipo, pi, w(po))^> 
^'^'^^ = -wi + w(Po) - 5po {ipi, {ipo + ptvo)^{'ipi + ptvi)) = 



where we have used (2.3). By the implicit function theorem the last 
equation has exactly one solution 

□ 

At Pq and at the corresponding value u* = u;(po), the family of even 
ground states which we have found above, bifurcates in two families, 
one formed by even ground states and the other by asymmetric ground 
states. In the context of the cubic NLS, |KKSWj proves that for u > u* 
the even ground states are unstable while the asymmetric ground states 
are orbitally stable. In the rest of §2 we double check that asymmetric 
ground states have the same behavior of |KKSW] for our quintic NLS. 

2.2. Asymmetric ground states for tu > tu*. Following |KKSWj we 
consider the branch of asymmetric ground states defined for u > u*. 
Set 
(2.7) 

F(po, pi, u) = -uopo + ujpo - {ipo, (Po^o + Pii^i + vipo, Pi, ^)y) , 

G(po, pi, oj) = -ui +UJ- — (ipi, (poV'o + PiV'i + viPo, Pu^)Y) ■ 

Pi 
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We know that F{pl,0,co*) = G{pq,0,co*) = 0. Hence, we apply the 
imphcit function theorem and prove the following 

Lemma 2.4. The Jacobian q\^^^^-^ has rank 2 at (pq,0,u;*). Corre- 
spondingly, by the implicit function theorem there are smooth functions 
Po(pi) o-nd uj{pi) such that 

^(Po(Pi),Pi,w(Pi)) = G'(po(pi),Pi,a;(pi)) = 0, 
which are defined for pi in a small neighborhood of and such that 



{21 



We have 



Po(Pi) = Po + 
u{pi) =u* + 



^"(0) = 20(p*) 



(2.9) 



Po(0) 



5 (v^o^v^?)-(v^o^v^^) 



pS5«V'?)-(^oM) 



5(^o',^?)-(V'oM) 
+ 0{{plf] 



+ o((pS)«)>o 



Proof. We have 



-Wo + w - 5(V^o, (Po^o + Pi^i 
+r/(po, Pi, w))^(V^o + dpoVipo, Pi, ^)))- 



For pi = and u; = ijj{pq) (see (2.3)), we get 

(2 \ {]\ = -uJo + uj- 5p^{ilJo, {^0 + PoVo)\^o + Po^i)) 

= -coo + co- 5p'M, 1) + 0{pl) = -ApM, 1) + Oipl). 

We have 

du^F = Po- 5(V^o, (Po^o + Pi^i + vipo, Pi,^)fd^v{po, Pu^)))- 
For pi = 0, we get 

(2.11) d^F = Po - 5p^(^o, (^0 + ptvo)''d^Vo)) = Po + O(p^). 



For pi = we already know from (2.5) that 

G = -uji + - 5po {^pi, iiJo + ptvo)^ii^i + PoVi)) ■ 

So, 

(2.12) a,„G = -20p^«,^i) + O(p^), a^G = l + 0(p«). 
Then the Jacobian matrix g^^pf^-^ at (po, 0, a;(po)) is 



(2.13) 



-4p^(^oM) + 0(pg) Po + 0(p9)\ 
-20pl{^lJ|,^l) + Oipl) l + Oipl)J 
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with inverse matrix 



20pg(V'o^ ^f) + o{pI) -^pti^l 1) + 0{pI) 

We see below that dp^F — dp^^G — ai pi — 0. This imphes immedi- 
ately a;'(0) = p(j(0) = 0. We have 

9p^F\^^^^ = -5(V'o, (PoV'o + PiV'i + ?7)''(V'i + 9pi77)) |^j=o 

= -5(^/'o, (poV'o + Po^)^(^i + Po^i)) = 0, 

where the last equality is due to the fact that ipi + PoTji is odd and the 
other factors are even. We have 

F\ 

- 5{^o, 4(poV'o + Pi^i + V?{^i + dp,vf 
+ (poV'o + PiV-i + vfdl,ri) |p^=o 

= -5(^0, 4(po^o + plmfi^i + Po^i)' + (poV'o + plr]of2plr]2) 
^-2Qpl{'il^l^l^l) + 0{pl). 

To show that dpiG\^^^^ — and compute ^piCj^^^g, we use the el- 
ementary fact that for two smooth functions f{x) = xg{x) we have 
^(0) = /'(O), ^'(0) = |/"(0) and /(O) = |/'"(0). Hence, calculating 

as above, 

^pMp,=o = lpi=0 (^1' (^0^0 + Pi^i + vf) = 

5 

- 2 (^i,4(poV'o + Po^o)'(^i + pIvi? + (Po^o + Po^o)'2p^r?2)> = 0, 

with the last equality due to the fact that we are integrating an odd 
function. Then, 

5 

- g^pi |^^=0 (^1' 4(^0^0 + P^^^ + + ^Pl'^)' 

+ (P0^0 + Pl^l+^)^f5pi^) 

= -^(V'l, 12(poV'o + Po^o)'(V'i + Po^i)' + (PoV'o + Po^o)^6p2r^3 

+ 12(poV'o + Po?7o)'(V'i + P^77i)2pg772) 
^-2Qpl{^l^l^l;{)+0{pl). 
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Hence we get, using M the matrix in (2.14) and the imphcit differen- 
tiation formula y' = —Qy^gx for g{x,y) = 0, and y" = —Qy^gxx when 
gx = 0, gy invertible, 

5 



□ 

2.3. The stability of the asymmetric ground states. We focus 
on 

0<^(pi) = Po(pi)^o + pi-ipi + v{po{pi),Pi, w(pi)), 
the asymmetric ground states. Following [KKSWj . we prove now that 
they are orbitally stable, that is for any such fixed u = uj{pi) and for 
any e > there is 5 > such that for uq G H^{M.) such that 

sup ||mo - e'^(f)^\\H^ < 6 



and for u{t) the solution of (1.1) we have that u{t) is globally defined 
and 

sup \\u{t) — e^'^cp^Wn'^ < e. 

Set 

and consider the pair of operators 
(2.15) 

L^iu) = -dl + V + u-bcl>t and L_ [u) = -dl + V + u - ct>l 

The proof of the orbital stability of asymmetric ground states is ob- 
tained in two steps. The second step is the following well known result, 
see [Wj, which in particular says that the e and the 5 in the definition 
of orbital stability can be here taken to be about of the same value. 
For the proof, see |Cu3j . 



Theorem 2.5. Given the hypotheses and conclusions of Lemma \2.6 
below, there 3eo > and Ao{(jj) > s.t. e G (0,eo) and 

||m(0,x) - < e 

imply for the corresponding solution 

inf{||M(t,x) - e*^0^(x)||jyi(M) : 7 G M} < Ao{uj)e. 
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The first step to prove the orbital stability of asymmetric ground 
states is the following 

Lemma 2.6. Consider the asymmetric ground states discussed in Sub- 
section\2^foruj > u* . Then, foru sufficiently close to u* the following 



two statements are true. 

(1) We have £q{uj) > 0. 

(2) For uj = w(pi) the set of eigenvalues of L^{uj) is given by 
ad{L+{uj)) = {-/io(pi),/ii(pi)}, where 

and /Ui(pi) > with fii{pi) ~ (Po)^Pi- 
Proof. We have 

^.M = ^(pg + P?) + ^(., = 2p.(l + 0(po)) 

+ :i^(Po°(^o, ^o) + PoP?(^i, ^i) + 2p^P?(^o, V2) + 0{pI)) 
"Pi 

= 2pi(l + 0(po)) + 0(p?). 



We have = a;"(0)pi + 0{pl) where u"{0) > by Q and ( fLSj ). 
Claim (1) follows from 

/ N dpi d 2pi(l + 0(po)) 2 , „ 

dj^^") = ^ = u"io)pi+oipi) = :m^' + > °- 

We consider the second claim. First of all, by ( |2.15[ ) and by the fact that 
is small, we know by general arguments that ad{L+{u)) is formed 
by two eigenvalues, the same number of (Trf(— 5^ + 1^ + 0;) (recall also 
that in dimension 1 the eigenvalues have always multiplicity 1). Since 

{^0, L+{co)^o) = (V^o, {-dl + V + u;)^o) - Hi^l <^t) 
= uj-ujo- 5(V^o, (poV^o + PiV^i + vY)- 



At = u*, p = p* and pi = 0, by (2.3) we obtain 

{ijjo, L+{u)^/jo) = - Wo - 5(^0' (Po^o + (Po)^'7o)'^ 

= -4(psr(v^o',i)+o((psr)<o. 

This means that L^[lj*) has at least one negative eigenvalue. By conti- 
nuity for 00 > u* close to u*, L^{u) has at least one negative eigenvalue, 
which we denote by —po{pi). We now start the discussion on pi{pi). 
We claim that 

(2.16) L4uj*)i^i + {pl)\{pl,uj*)) = 0. 
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u"{0)pi + 0((pi 



-(V'l + PoViipo, + 0{pi)) = -(j)^. 



Multiplying the equations by pi and letting pi \ we obtain (2.16). 
This implies that for small pi the eigenvalue problem L^{u))g = pg 
admits a small solution p = pi(pi) with pi(0) = 0. There is a unique 
solution L^{u)g = pg of the form g = ao'^o + i^i + ^, where ^ = 
Peg, p = pi and g are functions of pi. The expression L^{u)g = pg 
is equivalently expressed as follows considering a Lyapunov Schmidt 
reduction as in §5 in |KKSW] : 

(2.17) (^i,L+(u;)(aoV^o + V^i+6) =Ai, 

{~dl + V + io-p)^ = 5PXMo + ^1 + 0- 
We set = (Po)^'?^, where t]* = r]j{p^,u*). We have 

(2.18) [u-uo-p- 5(0^, ^o')) "0 = H<Pl Mi^i + 0)- 

By (2.16) we know a;o(Pi)|pi=o = 0. At pi = we have for u = u* 
(2.19) 



At pi = 0, from (2.17) and from ^'(0) = 0, see (2.8), we get 

20) ^''^ ^ 5Ryi-u*) [A<l>li^^ + e)2 + <l>t*a'oiJo + 

■ +p'iO)Rvi-co*)^* 

for Rv{z) = {-dl + V- z)-\ We have 

(2.21) p = u-uJi- 5ao{<Pt, V'oV'i) - H<Pt, Mi^i + 0)- 

Then, ( |2^ and ( |2^ imply p'(0) = 0. Indeed, at pi = we have 
uj'{0) = ao(0) = 0, hence we see 

p'(0) = -5a'o(0^,,^oV'i) - 20{<Pl,Mi^i+Q')- 



^id,,0 = -5p'(0)(</)^.,^ii?y(-u;*)e.). 



Since in (2.20) we have p'(0) = and by the fact that the resolvent 
Rv{z) preserves the spaces of even (resp. odd) functions, we conclude 
that dp-^C,\pi=o is even. We also have the estimate 

||9,,e|p.=o||L^<C^(Po)' + C(pS)Vo(0). 



This and (2.19) yield 



and so a'o(O) = -4(p*)-i^^^ 
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(^Jiy + '^(Po)- We have 
/i"(0) = u;"(0) -4Oa'o(O)(03(^, + 
(2.22) - 20(03^2^ 0^,7/>i(V'i + e*)) -60(02, (^, + ^j3^^^ 

Proceeding as above, |q;o(0)| < C{pl)~'^ and ||9^^^|p^=o||l°= 

< C. The 

terms in the second and third hnes of (2.22) are 0((po)^)- As a result, 
(2.23) 

/i (0) = 20(po) ^/ ,.4 ,.2\ 7^:^\ + 160(po) 



~6o{p*r{^plr,)+o{{plr)>o 



by (1.2H1.3). 



□ 



Lemma 2.7. yli pi = for the asymmetric branch we have the expan- 
sion 



(2.24) 



^ [Pi) 

4m(^o + 5(pS)%(pS,^*)) + (Po)'5.^o(pS,^*)- 
a; [U) 



Proof. By Lemma [TT] we have the following formula which follows from 

du(pu = -^dp,{po{pi)tpo + pi^i + r]{po{pi), pi,uj{pi))) 

dpi 



tIV"! + Por?i(po, w) + Po(^o + 5por?o) + uj'podu^Vo + 0{p^)). 

CO 

We the use Lemmas 12.41 

3. The discrete spectrum of the linearization 
Consider the operators 



□ 



(3.1) 



L_{io] 
-L+{uj) 



and J 



r 
-1 or 



Lemma 3.1. We have cfdii^u}*) = 0. For {} meaning span, we have 
(3.2) 



ker£^. = {ei{oo*), 62(00*)} with ei{oo) 







, 62(0;*) 



where we set P = ipi + (Po)^^?i- The generalized kernel is 
(3.3) N,{C.,) = {e,iu*):j = l,...A} 
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with 

esi^n = (o) ' = 

where (jyui* = L^{ijj*)a and (3 = L_(u;*)7 . We have that a{x) is an 
even function, while /3{x) and 7(x) are odd functions. 



Proof. The relation ( 3.2[ ) follows by the fact that is an eigenvalue of 



L_{u) and L+(a;*). The equations 0^* = L_^.{u*)a and /3 = L_{u*)'y 
admit solutions because {4>ui*,(3) = 0, since (f)^^* is even and (3 odd. In 



fact a coincides with the second line of formula (2.24). We conclude 



that Ng{Cuj*) contains the rhs of (3.3) and so dim Ng{Cui*) > 4. To see 



that they are equal it is enough to check that dim Ng{Cuj*) < 4. C^j* is 
a small perturbation of J(— <9^ + V + u*) which has 4 eigenvalues, all 
close to 0. This implies dim Ng^C^^*) < 4. □ 

Lemma 3.2. Let u > u* . Then ad{CiS) = {0, iA(a;), — iA(a;)} with 
X{u) > a simple eigenvalue. 

Proof. We know that d imiV q(£^) = 2, with Ng{C^) = {ei(cj), e2(w)} 



with 61(0;) as in Lemma 3.1 and with (e2(w))^ = {pid^cj)^, 0). Let D be 
a small disk containing the origin in its interior. We know that Ct^ — C^^* 
is continuous in u with values in the space of bounded operators from 
L^(]R) in itself with the uniform topology, and that C^* has exactly 4 
eigenvalues inside D (in fact just 1 with algebraic multiplicity 4) and 
that dD is in the resolvent set of C^j* . Then, 

ad(/:jnD = {0,iA(w),-iA(u;)} 

follows by the fact that ad{Cui) is symmetric with respect to the coor- 
dinate axes and has algebraic multiplicity 2 for C^., as we reminded 
above, and that the two eigenvalues cannot lie in M (since we know 
that the for uj > u* are orbitally stable). 

By standard arguments in perturbation theory, exploiting only 

\Mx)\<Ce-'^\-\ 

for a > and C > fixed and for both 11 = uj* , uj , it is possible to 
prove that 

ad{C^) n {C\D) 

is empty for oj close enough to u* . Specifically, and quite informally, 
elements of 

ad{C^) n {C\D) 

could originate by singularities of {C^^* — z)^^ on a second sheet of the 
Riemann surface where it is defined, or by the points ±iw* if was a 
resonance oi —d"^ + V . But the latter is excluded by hypothesis and 
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the former can be ruled out for uj close enough to u* . We skip the 
details: the analysis at the endpoints is similar to material in |Cu4j : 
the analysis of the eigenvalues coming from the second sheet can be 
derived from |CPVj . □ 



Lemma 3.3. Consider for cu > cu* the eigenvalue from Lemma 3.2 
with > 0. Let u > u* , then there exists a fixed C > such that 
A(w(pi)) > Cpip*. 

Proof. We recall that if C^U = i\U and = {u,v), then 

L^{uj)L^{ijj)u = \^u. 
Then, {u,(j)J) = and one can define / = which is s.t. 

(L_H)iL+M(L_M)5/ = AY 
Since one can proceed backwards, we have 

A — mm — 777: ^ 

mm —7^ mm ^--t^^ . 

{f,M=o \\f\\i2 {9,M=o \\g\\i2 

We prove now that 

Let 11/11^2 = 1 and = (j)^. Let L+{u)xj = (-)^^ViXi for j = 0, 1 
with llXjlU^ = 1- For a function g, let gj = {g,Xj) let gc = 
\\Pc{L^{uj))g\\L2, where (only for this proof) 

Peg ■= g - goxo - gixi- 

Then, 

/) > -/io/o' + /ii/i' + cof', =: F(/o, /i, /,). 

We will prove then that, subject to the constraints in the last two 
lines of (3.5) below, we have F > C2P1PI. Since F is continuous for the 
strong and weak topology in L^, there exists a constrained minimizer. 
This implies that, for a and b Lagrange multipliers, we have: 

2(10 - a)PJ = bP,^, 

-2po/o = 2a/o + 600, 

(3.5) 2pi/i = 2a/i + 60i, 

/o0O + /l01 + (Pc/,Pc0)=O. 
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For Pc0 and Pcf proportional to each other, the last equation in (3.5) is 
the same as /o0o + /i0i + /c0c = 0. If Pc4> and Pcf are not proportional, 
then 6 = and u = a. Then (w + /io)/o = implies fo = since 
u + fiQ > 0. Given {uj — fii)fo = 0, we have /i = since fj,i = 0{p\) 
while > > 0. Then, fc = 1 with F = u, which is clearly the 
maximum value, and not the minimum. Hence we can assume that Pc4> 
and Pcf are proportional. Then we minimize F under the constraint 



(3.6) 



0. 



/O0O + /l01 + /c</i 

Notice that the plane can be parametrized by /o = (piu + (pcV, fi 

-(pou, fc = -(pov- Then, 

(3.7) 

F(0iM + (j)cV, -(f)oU, -(poV) = -/io(0lM + (j)cV)^ + AH^o"^ + W0o^^ 
= (-/io0? + ^l4>l)u'^ + {-fJ'0(pl + ^^(Po)v^ - 2flo(j)i(f)cUV. 

This is a quadratic form in {u, v) with eigenvalues, x, the roots of 



(s - (-/io0? + fii4>l)){x - (-/io0c + ^4>l)) - 



(3.8) 

We have 

(3.9) 

and 

(3.10) 



= - (-Aio0? + /ii0o - A^o^c + ^0o) ^+ 



0. 



1 + 



(/ilW^o 



/iow0i - /ioAii0J0o- 
We claim that (3.10 )~ a;o(po)^pf. This and (3.9) imply that the polyno- 



mial in (3.8) has both roots positive, one about ujq{pqY and the other 
about {PqYp\- Then, the minimum of (3.7) for + f ^ = 1 is about 
{pIYpI- Since /q + /f + = 1 implies v/~+ ^ (Po)~^' follows that 
the minimum of F is > Clp^^pl for a fixed C. 

To show (3.10)~ a;o(Po)^Pi; '^^ observe that since L^{u) = L^{uj) — 
40fj, we have po < 4||0tj||^ < Cp^ for a fixed C. Then, the claim 
follows from 

I( |3l0| | > 0^(u;/ii0^ - Ca;p^P? - Cp^Pi^D 
2 

where we exploit 0i = 0(pi), 0o ~ Po; 0c = 0{pc) and pi 
p^p2. Then KM I >p^p2. □ 



> 00/^1 (^^00 



pIpI > 
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4. Set up for Theorem 11.21 and dispersion for the 

linearization 



Theorem 1.2 is a consequence of |Cu3] . Notice that since the hn- 
earization has just one pair of nonzero eigenvalues, the hamiltonian set 
up in |Cul] is unnecessary, and the theory in |Cu3t ICM] is adequate. 
We recall that due to absence of the endpoint Strichartz estimate in 
1 D, the theory requires some adequate surrogate. This for 1 D was 
provided by Mizumachi [M]. The theory in [M] though, is more compli- 
cated then necessary. The simplifications were provided in |Cu3t ICTj . 
Subsequent papers like jKPSl IPS ] return to more complicated approach 
[M]. So, even though Theorem 1.2 is a direct consequence of |Cu3j . we 
will take the opportunity to state the various steps of the proof, in 
order also to point out the points in jM] and in |PS] which can be 
simplified. 

Recall the ansatz 



(4.1) M(t,a;) = e^®W(</)^(t)(a;) + r(t,x)),e(t)= / co{s)ds + ^{t) 



Inserting the ansatz into the NLS (1.1) we get 

in = -r^x + Vr + uj{t)r - 30^i)r - 0^(i)f 
+ 7(^)0a;(t) - iw(t)<9c^0<^(t) + 7(t)r + 0(r 



2^ 



We set *i? = (r, r), *$ = {(p^, (p^) (using a different frame from the one 
in ^ and we rewrite the above equation as 

(4.2) iRt = H^R + asiR + (737$ - icod^^ + 0{R^) 

where: 

A / i\ /I 



^1 0) '^2= ^-i 0) '^3= \0 -1^ 
^^■^^ H^fi = (T^i-dl + V + uj),V^ = -3^3 [(Pt] + 2i0>2; 

We know that is an isolated eigenvalue of H^, dim Ng{H^) = 2. We 
have 

Since H* = a^H^a^, we have Ng{H*) = span{$(^, cr3(9<^$j^}. Let ^{00) 
be a real eigenfunction with eigenvalue X{u). Then we have 
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Notice that a^^) > since {(T^H^-, ■) is positive definite on N^{H*). 
For w G O, we have the iJ^-invariant Jordan block decomposition 

(4.4) L^{R,C^) = Ng{H^) (B {(B± N{H^T Huj))) (B Ll{H^), 

where Ll{H^) := {Ng{H*) © {®±N{H* t Huj))}^ ■ Correspondingly, 
we set 

(4.5) Rit) = zimcoit)) + zit)a,acoit)) + /(t), 

(4.6) R{t) e N^{H:^,^) and fit)eLl{H^^,)). 

There is a Taylor expansion at i? = of the nonlinearity 0{R^) in 
(4.2) with Rjn,n{^,x) and Am,n{^,x) real vectors and matrices rapidly 
decreasing in x: 0{R^) = 

2<m+?i<27V+l l<m+n<N 

Then, 
(4.7) 

ift = {H^(t) + / + a37$(w) - iw5^$(t) + {z\{w) - iz)i{w) 
- {z\{uj) + il)(Ti^(u;) + 0-37(2;^ + zaiS) - iLo{zd^^ + 2cri9^0 

+ z"'z^Rm,n{^)+ Y ■2™^"^m,n(c^)/+ 

2<m+n<2N+l l<m+n<N 

where by Oioc we mean that the there is a factor x(x) rapidly decaying 
to as |x| — )■ 00. Taking inner products of the equation with the 
generators of Ng{H*) and of ker(_fZ'* — A), we obtain modulation and 

discrete modes equations {q'{u)) := '^^^^"^^^ ): 
(4.8) 

iuq'iu) = {X, $) , 7g'(w) = {X, cTs^^^) , iz - \{uj)z = {X, cisO , 

2N+1 

X ■= a37(z^ + zaii) - '^^zd^i + zaid^S) + ^ z""z"i?„_„(w) 

m+n=2 

N 

+ (^37 + iud^Pc + J] ^"^^^^,,(0;))/ + 0(f) + 0,o,(|z2^+2|)_ 

m+n.=l 

We now go through the dispersive estimates. The proofs are in jCu3j . 
We call admissible a pair (p, q) s.t. 

(4.9) 2/p+l/g=l/2, p>4, g > 2. 
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Theorem 4.1 (Strichartz estimates). Fork = [0,2] there exist positive 
numbers C{uj,k,p) and C{u, k,pi,p2) upper semicontinuous in their 
arguments such that: 

(a) for any f G Ll{u)) and any admissible pair {p,q) with p > 4 we 
have 

(4.10) ||e-"^-/||^.^M<C||/||^.; 

(b) for any g{t,x) G ^(IR^) and any two admissible pairs {pj,qj) for 
j = 1,2 with pj > A we have 



(4.11) II / e-'(*-)^"P,(a;)^7(s,-)rfs||^.,^.„, < C\\g\ 
Jo * 



In the case k = 0, we can include also case p = A in (4.10) and pj = 4 
for any of j = 1,2 in (4.11). 

For the proof see |Cu3j . The case k > requires interpolation. 
The case /c = is hke the one for e~^*^^. Specifically, we can use 
dispersive estimates, see |KSj . and an appropriate version of the so 
called TT* argument. In particular, this yields the LfL'^ bound, which 
is not reached in [KS]. See |DMWj for how to extend such results 
to Schrodinger operators, H, formed by singular perturbations of the 
Laplacian with k < 1. 

Lemma 4.2. Fix r > 3/2. 

(1) There exists C = C{t,uj), upper semicontinuous in uj such that for 
any e 0, 

\\Rh^{X + ie)Pc{H^)u\\^2L^.-^ < C\\u\\l^. 

(2) For any u G L^'"^ the following limits exist: 

limRHSX ± ie)u = R% {X)u m C^{ae{H^), Ll'-^). 

(3) There exists C = C{t,u), upper semicontinuous in u such that 

ll^l(A)Pc(i^J|lB(L^.^,i^.-^)<C(A)-i. 

(4) Given any u G L^'"^ we have 

P,{H^)u = ^J (i?+JA) - Rj,^{X))udX. 

These are consequences of the fact that (Je{H^) does not contain 
eigenvalues and that ±aj are not resonances, and of the theory on 
plane waves and representation of the resolvent in |KSj . In fact, of 
a much simpler version than |KSj . due to the fact that H^^ is a small 
perturbation of cr^{—d1 + + a;). 
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Claim (a) of the following smoothing lemma is a consequence of 
Lemma 4.2 by [K], while (b) follows from (a) by duality. 



Lemma 4.3. For any k,T > 3/2, 3 C = C{T,k,u) upper semicontin- 
uous in u such that: 
(a) for anyfeSiR), 



Pc{H^)fh.Hi^-^<C\\f\U. 



(b) for any g{t,x) G ^(M^) 

e''''-P,{H^)g{t,-)dt 



Lemma 4.4. For any k,T> 3/2, 3 C = C{T,k,u) as above such that 
V g{t,x) G S{R^) 



f e-'^'-'^''-P,{H^)g{sr)ds 
Jo 



<C\\g\ 



2 rjk.T. 



LiH, 



Proof. To get this proof there is no need of Lemma 11 [M] or of the 
analogous result, Lemma 2 in §7, in |PS] . We just use Plancherel and 



r < 



Holder inequalities and (3) Lemma 4.2 
Jo * 



< \\RhJ^)Pc{Huj)xio,+oo) *a?(A, a;) 11^2^2.-. < 

< \\KA>')PciH^)\\B{Ll'\Ll'-^)\\X[0,+oo) *A?(A,X)||^2,. 

<II^H„(A)^c(^a;)|lL^(K,B(L2>-,L2> — ))||fi'|lL2i2,r < C\\g\\^2^2.r . 



/'2 



□ 



Lemma 4.5. k and r > 3/2 3 C = C{T,k,uj) as above such that V 
g{t,x) G S{R^) 



L^LlnLf{R,wi:''^) 



Jo 

Proof For g{t,x) G ^(M^) set 

r+oo 

Tg{t) = / e-"^'-'^''-P,{H^)g{s)ds. 
Jo 



Lemma 



4.3 



(b) implies / := e''^-Pc{uj)g{s)ds G L'^iR). Then 
Lemma [4. 5| is a direct consequence of |CKj . □ 



The following lemma can be proved in a way similar to Lemma B.l 
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Lemma 4.6. The following operators -P±(w) are well-defined: 



1 f^^ 

P+{uj)u = hm — hm / [Rh^ (A + le) - Rh^ (A - le)] udX, 

6^0+ ZTTl M-5>+oo 

1 r"" 

P {uj)u = lim — lim / [Rh^ (A + ie) - Rh^ (A - ie)] udX. 

e^0+ ZTTl A/-s>+oo 

For any M > and N > and for C = C{N, M, u) upper semicontin- 
uous in oj > u* , we have 

||(p+M - p_M - p,Ha3)/|U.,M < c\\f\\L^.-N. 

5. Normal form expansion 

Here we repeat the theory in |CMt ICu3] which is somewhat more 
elementary than |Culj . but still adequate in our setting. 

We consider iV G N such that for any t > 0, for pi{t) and for the 
corresponding u(t) := uj{pi(t)) and for A(t) = \{uj{pi{t))) we have 

(5.1) NX{t) < uj{t) <{N + l)A(t). 

Notice that '^^^l^p^'l^ is for pi > a continuous and strictly increasing 
function, equal to at pi = 0. This means that for pi > small, it 
has no values in N. By continuity and orbital stability, we can then 
assume (5.1) for all t. 



5.1. Changes of variables on /. For the of (5.1) we consider 
k = 1,2, ...N and set f = fk ior k = 1. The other are defined below. 
In the ODE's there will be error terms of the form 

EoDEik) = 0(|zp^+2) + 0(z^+Vfc) + 0{fl) + 0(|/,n. 

In the PDE's there will be error terms of the form 

EpDEik) = 0z„,(i^r+2) ^ Oloeizfk) + oufi) + om'fk). 

In the right hand sides of the equations (2.3-4) we substitute 7 and Co 
using the modulation equations. We repeat the procedure a sufficient 

number of times until we can write for = 1, /i = / and q'{ijj) = 
(5.2) 

2N+1 N 

iuq'iu) = ( J2 ^'"^"AWJo;) + Yl z^-z-Ati{uj)f, 

m+n=2 m+n=l 

+ EoDE{k),^{uj)) 
iz — \z = { same as above , 0-3,^(0;)) 

idtfk = (H^ + ^37) fk + EpDE{k) + ■^"^"^L'U^^), 

k+l<m+n<N+l 
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with Am]n, Rrn,n and Am,n(w,a;) real exponentially decreasing to for 
|a;| — 7- cxD and continuous in {oj.x). Exploiting \{m — n)A(a;)| < u for 
m + < A^, m > 0, n > 0, we define inductively fk with A; < by 

fk = z"'z''^mA^) + fk-i for "^mA^) 

m+n=k 

Notice that if RmA\'^i ^) is real exponentially decreasing to for |x| — )■ 
oo, the same is true for "i/mA^) \{m — n)A(a;)| < u. By induction 
fk solves the above equation with the above notifications. 

We are now ready to state the result which directly implies Theorem 
Ol 



Theorem 5.1. Assume (H1)-(H5). Let u he a solution of ), U 



\u,u), and let "^rnA^) '^^ above. Then if to in Theorem 1.2 is 
sufficiently small, there exist -functions u{t) and 0{t), a constant 
> (jj* such that we have supi>o — ujq\ = 0{e), lim(_j.+oo ^(t) = 
u;+ and we can write 

U{t,x) =e^^W'^^ + z{t)i{At)) +W)^ii{At))) 

+ j2 ^,^AAt)HtrWf + e''^'^'^^f^{t,x), 

2<m+n<N 

with + \\fN{t,x)\\^^^i^^r 1,10 ^4^^ < Ce. 

Furthermore, there exists /+ G i/^(]R, C^) such that 

= 0. 



(5.3) lim e^''W"VA'(i) - e'*^'"'/ 

>+oo 



Obviously the scattering result (5.3) holds also for t — )■ — oo. We 



do not prove this lemma explicitly, but we recall Lemma 4.3 in |Cu3] 
which states: 

Lemma 5.2. There are fixed constants Cq and Ci and eo > such that 
for any < e < eo if we have 

(5-4) \\z\\^p\2 < 2Coe & \\fN\\LocfjinLlwt'°nLiL^n^Ht-'^'^^'^'^^ 
then we obtain the improved inequalities 

(5-5) ll/A^llL,°°HinL5iyi'i°nL4L^nL2/fi>-2 < C'le, 

(5.6) ||2||^+V2 < Coe. 
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We sketch only the main steps of the proof. First of all we rewrite 
the equation for /jy. Set c<;(0) := c<;(pi(0)) and write 

idth - {H.io) + as (7 + - u;(0))) Jn = J] 

m+n=N+l 

+ Epde{N) where Epde{N) := Epde{N) + (V^ - V^(o))/jv. 
It is easy to see that (5.5) for f^, which is the same of Pc{u!)fN, 



or for Pc{uj{0))fN, are equivalent. This because Pc{oo) — Pc{oj{0)) = 
Pfi{u}{0)) — Pd{uj) is a small and smoothing operator. We then write for 

ip(t) = ip = 'J + u — u{0) 

idtPMo))fN - {H^m + V {P+{^m - P-{^m)) PMo))fN = 

J2 ^'"^"Pc(c.(0)X)(u;)+ 

m+n=N+l 

P,{ujmEpDE{N) + ip (P+(a;(0)) - P-(a;(0)) - PMO)M In- 
It turns out that the term in the second line is the main one of the rhs 



and that the norms of /at in (5.5) can be controlled by ||/Af(0)||j:/i + 



lkll'!Sv+2- In particular for Epde{N) we refer to Lemma 4.5 |Cu3j . 
The very last term in the equation is controlled using the fact that ip 



is small, Lemmas 4.4, 4.5 and 4.6 



5.2. A further change of variable in /. In the argument there is 
need for a decomposition of /at, namely setting 

(5.T) fN = - ^™^"^t(o)(("^-^)^)^^(^(0)X"H+^' 

m+n=N +1 

where if |A| < a;(0), we set R^j ^qj(A) = RH^^^y^{^)■ The function g 
satisfies an equation of the form 

idtP,{ujmg = {H^^o) + V^{P+{iom - P-HO)))) PMO))g+ 
■sr^ i \N+i\ 



+ ^ 0(e|zr^+^)i?+^^^,^ (±(iV + l)A(a;(0)))i?± + P^^umEpoE^N) 
± 



R+ := P);i,o, R. := and EpBs(iV) := Epde{N)+OU^z''+^). 



Lemma 5.3. Assume the hypotheses of Lemma\5.2 Then, there exists 



a fixed Co = C{u{0)) such that for a fixed S sufficiently large 
(5.8) \\gh,ri,-s<Coe + 0{e'). 



See Lemma 4.6 |Cu3j . 
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5.3. Change of u and z. Consider now the equations of uj and z in 



(5.2). Then, we have the following 
Lemma 5.4. There is a change of variables 

u = u + q{u,z,z)+ ^ 2;"'f"(/7v, A„„(w)), 

(5 9) l<m+n<N 

C = Z+P{UJ,Z,Z)+ ^'"^"(/^.i^-nM), 
l<m+n<N 

with p{u, z, z) = J2Pm,ni^)z"^z"- andq{z,z) = J2'lm,n{^)z"^z"- polyno- 
mials in {z,z) with real coefficients and 0(|2;p) near 0, such that we 
get for amioj) real 

itj = {Epde{N),<!>) 
(5.10) iC-AMC= Yl (^m{com'"'C+{EoDE{N),asO 

l<m<N 

Proof. The proof is elementary and goes as follows, see |CMj for details. 
We consider recursively for i = 0, ...,2N with zq = z equations 

l<l<i l+2<m+n<2N+l 

+ Y zTz-dA['U^)^fN)+EoDE{t) 
e+l<m+n<N 

with aiilniuj) e R and A^ln(w) G Suppose this holds for 

i < 2N. Then set 

m+n=e+2 m+n=£+l 
W ( \ 

■■= , "T^^^^, for m ^ n + 1, /^S.Jc) = 0, 
(m — 1 — n)A 

Bl^liu) = -i?H*((m - 1 - n)X)Aliliuj) for i < N, 

Bl^liuj) = -RH*{{m - 1 - n)X)Ai^]^{u) for {m,n) ^ (0,iV), 

where i?Q^(u;) = otherwise and Bm]n{uj) = for i > N. This yields 
for ( = Z2N the desired result. 



We substitute in the equation for uj in (5.2) (for k = N) z with 



( inverting the second equation in (5.9). We consider recursively for 
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i = 0, 2N + 1 with Qq = u equations 

e+2<m+n<2N+l e+l<m+n<N 
+ EoDEii) 

with 7£„(w) G M and r5n(w) G 5(M^M2^. Suppose this holds for 
i <2N + 1. Then set 

m+n=£+2 m+n=e+l 

AWJcu) = - n)X)Tliliu) for £ < N, 

Am,n{uj) = for i > N. This yields for oo = ^2N+i the desired result. 

□ 



By (jKij) we have = O(e^). 

Remark 5.5. Setting cD(t) = a;(0), /^r = and considering the equation 

iC-A(a;)C= E «-HICr"C 



l<m<Af 



yields a finite dimensional approximation of the NLS. We do not check 
here the time span when the solutions of this approximation are good 
approximations of solutions of the full NLS. Nonetheless we recall that 
in the literature, see for example |BP2t IGSj . are displayed solutions 

K(o)| 



of (5.10) s.t. approximately \Cit)\ ~ q-, with this 

^ — ^ (|C(o)|2^ivrt + i)27v' 

approximation valid for t <^ |C(0)|~^^. 



6. The Fermi golden rule 



Substituting in the equation for C the variable /tv with (5.7) to get 
iC-A(a;)C= E (^MlCrC +{EoDEiN),crsO- 



l<m<N 

h 
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With a^, <^ and R^^},^, real. Set 

T{u,um = Im ((45Mi?+^^J(iV + l)A(a;(0)) 

= 7r(<^M5(i7^(o) - (iV + l)A(a;(0)))Pe(a;(0))i?!v'2i,oM^3eM)- 

Now we assume the following: 

(H5) There is a fixed constant F > such that |r(a;,ci;)| > F. 

It is then easy to see, for example Corollary 4.7 |Cu3] . that in fact 
r(ci;,a;) > F, but we will not use this here. By continuity, we can 
assume |F(co', Ci;(0))| > F/2. Then, we write 

1^ = -F(c.,a;(0))|Cp^+^ + Im((<iM/^+i,a3eM)r^^ 
+ lm{{EoDE{N),a3a^))C). 



For Aq an upper bound of the constants Ao{u) of Theorem 2.5, we get 



Then we can pick Cq = 2(Aq + 2c{u{0) + 1))/F in Lemma 5.2 and this 
proves that (5.4) implies (5.6). Furthermore ((t) — )■ since ^Cit) = 



0(e), again see " fCuSj for more in depth discussion 

Appendix A. Finite Dimensional Dynamics 
We plug the ansatz 

u{x, t) = co{t)ipo + Cl^pl + R{x, t) 



into (1.1), where 

for j = 0, 1. As a result, we have the equation 

icoi>o + iciijji + iRt{x, t) = -uoCoipo - UiCiipi + HR 

-\coipo + ciipi + Rl'^icoipo + ci^i + R). 
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The nonhnear contribution is then given by 

{coipo+cM%CotPo + CitPif + 0{R) = 

+OiR) = 

{clcl)ip^+{3clclci + 2clcoCi)ipQ'4)i + (ScoCqCi + 6clcoCiCi 

+clcl)^l^l + {clcl + QcoCocjci + SclcicD^lipl 

+ (2coC?Ci + ScoclcD^Po^t + (clcD^Pl + 0{R). 

Projecting onto ipQ and ipi respectively and for now ignoring com- 
ponents with dependence upon R (see |MWj ) . we arrive at the finite 
dimensional Hamiltonian system of equations given 

ipo(V^o,^o) = -WoPo(V^o,^o) -PoPo(^O'V'o) 

-(SpoPoPi + 6poPoPiPi - PoPi)(^0' ^i) 
-(2poP?Pi + 3poP?P?)(^o>V^t), 

^Pi(^i,V'i) = -wiPi(V^i,V^i) -p?p?(V^i,V'i) 

-(SPoPoPi + 2poPoPi)(^0'^?) 

-(6poPoP?Pi + SpoPiP? + PoP?)(V'o5 V't) 

and the corresponding conjugate equations. Note, mass is conserved in 
this finite dimensional system, hence we have 

\po\'' + \pi? = N 

for all t. 

Plugging in alternative coordinates designed to give rise to a simple 
classification of the finite dimensional dynamics, we set 

Po(t) = A{t)e'^'^ 

and 

pi(t) = (a(t) + z/3(t))e^^W. 
As a result, we have 
iA-Ae = -uqA- A^ -3A^{a + i/3f -6A%a^ + 13^) 

-2A{a + i(3f{a' + (3^) - A\a - i(3f - ?,A{a' + (3^f 

and 

i{a + 0) - {a + il3)9 = -uJi{a + iP) - 3A^{a^ + (3^){a - 0) 

-A\a + ipf -2A\a-tP) 
-6A\a^ + P^){a + i(3) - ?>A\a + i(3) 
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setting for simplicity = 1 for all j = 0,1,2,3. This will 

simply rescale the dynamical system and not impact the general shape 
of the phase diagram. 
In the end, we have 

a = (a;o-^i+4v4V + 2(a2 + /32)2 + 2(a2 + /32)(a2-/32))/5, 
$ = -(wo -c^i-4A^- 4^2/32 + AV + 2(a2 + /32)2 

A = -AA{A^ + {a^ + P^))aP, 

e = wo + ^' + 10AV + 2A2/32 + 2(«2^/32)(a2_/32) 
where 

N = A^ + a'^ + /3l 

Using the mass conservation, we can write a closed system for [a, 
From the equation for /3, it is clear that in this rescaled dynamical 
system, we have 



FD 



LUO — CUI 



4 

We observe in Figure [T| several phase diagrams for varying values of 
A^, which point out the existence of periodic solutions above, near and 
below the bifurcation point. It is our goal in this section purely to 
give further evidence that the quintic NLS with double well potential 
presents similar dynamics to that of the cubic NLS with double well 
potential. For a dynamics approach to classifying these solutions and 
studying their stability properties, we refer to the finite dimensional 
results in |MWj for techniques which directly apply to reducible Hamil- 
tonian systems of this type, particularly for the proof of existence of 
periodic orbits and the resulting Floquet stability analysis. However, 
as the intent of this note is to prove asymptotic stability, we do not 
explore this topic further here. 

Appendix B. Numerical Verification of Hypotheses 
For a potential well 

V{x) = (j)„{x - L) + (j)^{x + L) 

with 
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Figure 1. Symmetric state for - N^^^ < corre- 
sponds to an equilibrium elliptic point {a, (3) = (0,0) 
and asymmetric states for — A*"^-^ > correlate to 
equilibrium points at {a, (3) = {±acr,0). Plotted phase 
portrait corresponds to parameter values ojq — ui = .1 
and = .1, .2, .5 respectively. 



where a = .001, L = 7.5. We discretize using a finite element method as 
in |MW] on a finer and finer set of grids that are concentrated near the 
peaks of the delta functions, we compute using the standard eigenvalue 



and eigenfunction solvers from Matlab the values from (1.2), (1.3) from 
hypothesis (-^4) are computed as 



and 



5(7/.o'>V^?)- (^0,1) ~ -3305 



20 .o , +160 ^,°'/,\^ - 60(^, 



5(^o',^?)-(^oM) 



9.9143 



respectively, showing that computationally at least our hypotheses are 
valid for a particularly interesting symmetric potential. 
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